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ABSTRACT 


Numerical  solutions  of  the  compressible  Navier-Stokes  equations  are  "presented* 
for  a  laminar  horseshoe  vortex  flow  created  by  the  interaction  of  a  boundary  layer 
on  a  flat  surface  and  an  elliptical  strut  leading  edge  mounted  normal  to  the  flat 
surface.  The  computational  approach  utilizes  "zone  embedding*"?^ surface-oriented 
elliptic-cylindrical  coordinates,  interactive  boundary  conditions,  and  a  consistently- 
split  _linearized  Wock  implicit  (LBI)  scheme  developed  by  the  authors.  Mesh  resolu¬ 
tion  tests  are  performed,  and  the  horseshoe  vortex  flow  is  discussed. 

\ 

INTRODUCTION  \ 


The  present  study  considers  three-dimensional  leading-edge  horseshoe  vortex  flow 
representing  the  interaction  of  a  two-dimensional  wall  boundary  layer  approaching  an 
elliptical  strut  leading  edge  mounted  normal  to  a  flat  surface.  An  example  of  this 
type  of  flow  occurs  near  the  junction  of  an  airfoil  or  wing  and  its  supporting  sur¬ 
face.  Another  example  is  present  in  axial  compressors  and  turbines,  where  boundary 
layers  which  develop  on  the  annular  surfaces  of  the  axial  flow  passage  encounter  rows 
of  stationary  and  rotating  blades.  The  flow  is  of  interest  in  connection  with  flow 
degradation,  its  tendency  to  cause  high  local  heat  transfer  rates,  and  its  role  as 
the  origin  of  corner  flows.  Numerous  experimental  flow  visualization  studies  (e.g., 
[1-3])  have  established  that  the  flow  consists  of  a  three-dimensional  boundary  layer 
separation  in  front  of  the  obstruction  followed  by  a  vortex  flow  which  wraps  around 
the  obstruction.  Little  is  available  in  the  way  of  detailed  flow  measurements,  how¬ 
ever,  particularly  downstream  of  separation.  Previous  analytical  studies  have  con¬ 
sidered  the  three-dimensional  boundary  layer  flow  upstream  of  separation  (e.g.,  [4]) 
and  have  used  rotational  inviscid  flow  theory  to  estimate  secondary  flows  [5].  The 
horseshoe  vortex  flow  is  treated  here  by  numerical  solution  of  the  compressible 
Navier-Stokes  equations  for  laminar  flow  at  moderate  Reynolds  number  Re  and  low  Mach 
number  M.  The  study  has  as  its  goals  an  improved  understanding  of  the  horseshoe 
vortex  region  and  formation  of  corner  flows,  and  also  the  development  of  flow  pre¬ 
diction  techniques. 


ZONE  EMBEDDING  AND  INTERACTIVE  BOUNDARY  CONDITIONS 


The  compressible  Navier-Stokes  equations  in  general  orthogonal  coordinates  are 
solved  using  analytical  coordinate  data  for  an  elliptic-cylindrical  coordinate  sys¬ 
tem  which  fits  all  solid  surfaces  within  the  computational  domain  but  is  not  aligned 
with  the  direction  of  the  free  stream  flow.  In  selecting  the  computational  domain, 
a  "zone  embedding"  approach  is  adopted  whereby  attention  is  focused  on  a  subregion 
of  the  overall  flow  field  in  the  immediate  vicinity  of  the  leading  edge  horseshoe 
vortex  flow.  A  perspective  view  of  the  geometry,  coordinate  system  (5,n,z)  and  a 
representative  computational  grid  is  shown  in  Fig.  1.  The  elliptic  coordinates  5 , n 
are  related  to  Cartesian  coordinates  x,y  by  x=cosh  5  cos  n,  y*sinh  5  sin  n,  and  the 
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(conformal)  metric  scale  factor  h  is  given  by  h  =cosh  5-cos  n.  The  computational 
domain  is  defined  by  0.2£5sl.6,  *r/ 2£ni3*r/2,  0<z<l. 

Symmetry  conditions  are  applied  at  z-1,  so  that  the  flow  represented  is  that 
past  an  elliptical  cylinder  (of  major  semi-axis  a«1.02  and  minor  semi-axis  b*0.2) 
between  parallel  flat  plates  with  spacing  H"2.  Velocity  boundary  conditions  at  no¬ 
slip  and  symmetry  surfaces  are  straightforward  and  self-explanatory.  The  remaining 
condition  applied  at  these  boundaries  is  3p/3n-0,  where  p  is  pressure  and  n  denotes 
the  normal  coordinate  direction.  The  condition  3p/3n»0  at  a  no-slip  surface  is  cor- 


rect  to  order  Re  for  viscous  flow  at  high  Reynolds  number. 

The  treatment  of  inflow  and  out¬ 
flow  conditions  is  the  principal 
obstacle  to  be  overcome  within  the 

present  zone  embedding  approach.  At  Z 

curved  boundaries  located  within  the 
free  stream  region,  interactive 
boundary  conditions  are  derived  from 
an  assumed  flow  structure  and  physi¬ 
cal  approximations  which  permit  in¬ 
flow  of  a  boundary  layer  and  inviscid 
free  stream,  and  which  permit  outflow 
in  the  presence  of  shear  layers,  cor¬ 
ner  flow,  streamwise  vorticity  and  a 
nonuniform  free  stream.  The  initial 
conditions  and  interactive  boundary 
conditions  are  devised  from  the  in¬ 
compressible  potential  flow  velocity  Fig.  1  -  Geometry,  Coordinate  System,  and 
Uj(C,ri)  about  the  ellipse,  together  Representative  Grid 

with  two-dimensional  estimates  of  the 
boundary  layer  thicknesses  on  the  end- 
wall  flat  plate,  6^(x),  and  on  the 
ellipse,  and  finally  from  an 

estimate  of  the  blockage  correction  factor  B(x)  for  the  core  flow  velocity  due  to  the 
endwall  boundary  layer  growth.  The  complex  potential  W  for  flow  past  an  ellipse  at 
incidence  a  with  circulation  k  is  given  by  Milne-Thomson  [6]  as  W=Ur(a+b)  cosh 

(£-£  -ia),  where  U  is  a  reference  free  stream  velocity,  £=£  is  the  ellipse  surface, 

°  r  2222° 
and  ?  is  defined  by  x+iy=c  cosh  £,_where  i  =~land  c  *a  -b  .  This  defines  the  incom¬ 
pressible  potential  flow  velocity  UT.  The  boundary  layer  thicknesses  6.  and  6„  are 

L  1/2  1  1 

approximated  by  the  Blasius  value  5(s/Re)  where  s  is  distance  from  the  flat  plate 

or  ellipse  leading  edges,  as  appropriate.  The  blockage  factor  is  given  by  B(x)= 

[H/2-6*(x) ]-l.  Finally,  boundary  layer  velocity  profile  shapes  f^Cz/5^)  and 

f2(C><$2)>  are  de^ined  from  vonKarman-Pohlhausen  polynomial  profiles.  The 

initial  velocity  vector  0  at  t=o  is  defined  by 


U(£ ,tj,  2  )  =  Ujte,??)  B(x)  f,  [  Z/8|(X)]  f2[(  C-£0)/  M7?)] 


(1) 

At  the  inflow  boundary  a  "two-layer"  boundary  condition  is  devised  such 

that  stagnation  pressure  pQ  is  fixed  at  the  free  stream  reference  value  (PQ)r  in  the 

core  flow  region  (z>6^),  and  the  Cartesian  velocity  u  is  set  by  u“ug(n,t)f1(z/6^)  for 

z<6, .  Here,  u  is  the  local  free  stream  velocity  consistent  with  p  and  the  local 

wall  static  pressure  (assumed  constant  across  the  shear  layer),  which  is  determined 
as  part  of  the  solution  and  updated  after  each  time  step.  The  remaining  inflow  con¬ 
ditions  are  v-aU^Cz/5^ ,  32w/3n2«o,  and  32cp/3n2-g(n,z)  where  g  is  the  distribution 

of  this  quantity  at  t“o  with  cp  defined  as  (l-B^U^Uj),  its  value  from  the  potential 

flow  corrected  for  estimated  blockage.  The  velocities  v  and  w  are  the  Cartesian 
velocity  components  normal  to  the  plane  of  the  ellipse  chord  and  normal  to  the  flat 
plate,  respectively.  The  angle  a  is  the  approach  flow  incidence  angle  relative  to 
the  chord  of  the  ellipse.  For  outflow  conditions,  second  normal  derivatives  of  each 
Cartesian  velocity  component  are  set  to  zero  and  cp  is  imposed  and  updated  after  each 


time  step  from  an  interaction  model  relating  the  computed  outflow  velocities  and  the 
a  priori  potential  flow  solution.  The  imposed  cp  is  obtained  by  integrating  an  in- 

■} 

viscid  normal  momentum  equation  3p/3n=-q“/Rj  along  the  outflow  symmetry  lines 

(z=l,n“f/2,37r/2)  beginning  with  a  fixed  wall  pressure.  Here,  q  is  the  local  computed 

2_-  - 

tlow  speed  (q  =U-U)  and  is  the  local  potential  flow  streamline  radius  of  curvature 
The  resulting  static  pressure  distribution  is  imposed  at  the  outflow  surface  assuming 
no  variation  normal  to  the  endwall. 

The  foregoing  interactive  inflow-outflow  boundary  conditions  are  designed  to 
permit  the  mass  flux  through  the  computational  domain  to  adjust  to  both  the  imposed 
downstream  static  pressure  and  to  viscous  losses  present  in  the  flow,  while  maintain¬ 
ing  a  specified  flow  structure  based  on  physical  assumptions  consistent  with  the  flow 
problem  under  consideration.  Various  refinements  in  the  interactive  boundary  condi¬ 
tions  are  possible,  such  as  including  the  effect  of  local  pressure  gradient  on 
boundary  layer  growth  and  profile  shape. 

DIFFERENCING  PROCEDURES 

The  differencing  procedures  used  are  a  straightforward  adaptation  of  those  used 
by  the  authors  [7]  in  Cartesian  coordinates  for  flow  in  a  straight  duct.  The  com¬ 
pressible  time-dependent  Navier-Stokes  equations  are  written  in  general  orthogonal 
coordinates,  and  for  economy  the  stagnation  enthalpy  is  assumed  constant.  The  defini¬ 
tion  of  stagnation  enthalpy  and  the  equation  of  state  for  a  perfect  gas  can  then  be 
used  to  eliminate  pressure  and  temperature  as  dependent  variables,  and  solution  of 
the  energy  equation  is  unnecessary.  The  continuity  and  three  momentum  equations  are 
solved  with  density  and  the  C,n  and  2  velocity  components  as  dependent  variables. 
Three-point  central  differences  were  used  for  spatial  derivatives,  and  second-order 
artificial  dissipation  terms  are  added  as  in  [7]  to.  prevent  spatial  oscillations  at 
high  cell  Reynolds  number.  This  treatment  lowers  the  formal  accuracy  to  first  order 
but  does  not  seriously  degrade  accuracy  in  representing  viscous  terms  in  thin  shear 
layers.  Analytical  coordinate  transformations  were  used  to  redistribute  grid  points 
and  thus  improve  resolution  in  shear  layers  and  near  the  leading  edge.  Derivatives 
of  geometric  data  were  determined  analytically  for  use  in  the  difference  equations. 

SPLIT  LBI  ALGORITHM 

The  numerical  algorithm  used  is  the  consistently-split  "linearized  block  im¬ 
plicit"  (LBI)  scheme  developed  by  the  authors  [7,  8]  for  systematic  use  in  solving 
systems  of  nonlinear  parabolic-hyperbolic  partial  differential  equations  (PDE's). 

To  illustrate  the  algorithm,  let 


(<£"*'  -  <£")/At  =  pD  (<£n  +  l)  +  (l-/3)D(*n) 


(2) 

approximate  a  system  of  time-dependent  nonlinear  PDE's  (centered  about  tn+6At)  for 
the  vector  J  of  dependent  variables,  where  D  is  a  multidimensional  vector  spatial 
differential  operator,  and  t  is  a  discretized  time  variable  such  that  At=tn+l-tn.  A 
local  time  linearization  (Taylor  expansion  about  $n)  is  introduced,  and  this  serves 
to  define  a  linear  differential  operator  L  such  that 


D(  ^n*')  =  D(<£n)  +  Ln  (<£"*'-  <£n)  +  0(  At2) 


Eq.  (2)  can  thus  be  written  as  the  linear  system 


(3) 


(I  -/?AtL n)(<£n>'-<£n)  =  A t  D ( <pn) 


(4) 


The  multidimensional  operator  L  is  divided  into  three  "one-dimensional"  sub-operators 
L=L1+L,+L3  (associated  here  with  the  three  coordinate  directions),  and  Eq.  (4)  is 

split  as  in  the  scalar  development  of  Douglas  &  Gunn  [9]  and  is  written  as 

(r -/3AfL,n}(<£*  -  <fin)  =  Ato(<£n)  (5a) 

(I  -jS  At  L2n  ){^>**  -  <fP)  =  4>*-4>n  (5b) 

(I  -/SAtL3n)(<^***  -  4>n)--  <£**-<£"  (5c) 

<£"*'  =  <£***  +  0(Af3)  (5d) 

If  spatial  derivatives  appearing  in  L  are  replaced  by  three-point  difference  formulas, 
then  each  step  in  Eqs.  (5a-c)  can  be  solved  by  a  block-tridaigonal  "inversion". 
Eliminating  the  intermediate  steps  in  Eqs.  (5a-d)  results  in 


(I  -/3At  L,n)  ( I  -  /3  At  |_2n  )  ( I  -/3AtL3n)(^n+l-^n)  =  At  D(0n) 

(6) 

3 

which  approximates  Eq.  (4)  to  order  At  .  Complete  derivations  are  given  by  the 
authors  in  [7,  8].  It  is  noted  that  Beam  &  Warming  [10]  have  reformulated  this  al¬ 
gorithm  as  a  widely-used  "delta  form"  approximate  factorization  scheme  whose  two- 
level  form  is  identical  to  Eq.  (5). 

MESH  REFINEMENT  IN  TWO  AND  THREE  DIMENSIONS 

Solutions  were  computed  for  three-dimensional  horseshoe  vortex  flow  past  an  el¬ 
liptical  leading  edge  for  free  stream  angle  of  incidence  a  of  zero  and  five  degrees. 
These  solutions  have  the  following  flow  parameters:  5:1  ellipse,  chordal  Reynolds 
number  of  400,  Mach  number  0.2,  and  flat  plate  leading  edge  located  1.25  chords  up¬ 
stream  of  the  leading  edge  of  the  ellipse.  For  reasons  of  economy,  these  solutions 
were  computed  using  a  coarse  computational  grid,  the  finest  being  14x28x14  ($,n,z) 
for  i=5°.  However,  the  geometry,  flow  conditions,  and  grid  distributions  were  care¬ 
fully  adjusted  to  provide  the  best  resolution  possible  within  this  constraint.  Most 
length  scales  expected  to  be  present  in  this  type  of  flow  receive  at  least  modest 
resolution.  Specifically,  care  was  taken  to  provide  resolution  of  boundary  layers  on 
the  endwall  and  ellipse,  of  the  Heimenz  layer  at  the  leading  edge  stagnation  point, 
and  of  the  flow  region  within  one  leading  edge  radius  of  both  the  leading  edge  and 
the  endwall.  The  only  potentially-relevant  unresolved  length  scale  is  an  0(Re~3/8) 
streamwise  distance  predicted  by  two-dimensional  "triple  deck"  theory  near  a  separa¬ 
tion  point. 

The  truncation  error  associated  with  the  present  flow  and  computational  grid 
was  examined  by  mesh  refinement  and  by  comparison  with  an  incompressible  solution  of 
Lugt  and  Ohring  [11]  for  two-dimensional  flow  past  a  10:1  ellipse  at  zero  incidence 
and  Re=200.  The  latter  comparison  for  streamwise  velocity  at  the  present  outflow 
boundary  (p=t/2)  is  shown  in  Fig.  2.  The  agreement  is  quite  reassuring  in  light  of 
the  respective  differencing  procedures  and  meshes.  Good  agreement  was  also  obtained 
for  surface  pressure  distributions.  Additional  two-dimensional  solutions  for  the 
forward  half  of  a  5:1  ellipse  (a**0°,Re=200,q=0.2)  were  computed  assuming  symmetry 
about  the  chord  line  n*ir  and  using  both  14x14  and  28x28  grids.  The  streamwise 
velocity  U  along  the  stagnation  streamline  is  a  sensitive  indicator  of  mesh  depend¬ 
ence  and  is  shown  on  a  logarithmic  scale  in  Fig.  3.  The  14x14  mesh  is  identical  to 
the  14x28  mesh  used  for  £»n  in  the  three-dimensional  solutions  without  symmetry  about 
n*ir,  and  the  small  amount  of  mesh  dependence  in  Fig.  3  indicates  that  this  mesh 
should  be  adequate,  at  least  away  from  the  endwall.  Finally,  resolution  normal  to 
the  endwall  (z  direction)  was  examined  for  three-dimensional  flow  at  5°  incidence 
using  both  14x28x10  and  14x28x14  grids.  The  velocity  profile  at  x»-1.06,v*0  near 
separation  is  the  most  sensitive  indicator  and  is  shown  in  Fig.  4.  The  finer  mesh 
improves  resolution  considerably  very  near  the  endwall  flow  reversal,  but  makes 
little  difference  elsewhere. 
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Fig.  2  -  Velocity  Profiles  at  Outflow 
Boundary  ( n=  n / 2 )  for  10:1 
Ellipse",  a=0°,  Re=200. 
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The  foregoing  results  provide  a 
limited  but  informative  assessment  of 
truncation  error  associated  with  the 
present  calculations.  The  solutions  may 
be  characterized  as  having  both  qualita¬ 
tive*  significance  and  reasonable  quanti¬ 
tative  accuracy.  Any  potential  concern 
that  failure  to  provide  high  resolution 
near  separation  would  cause  serious  er¬ 
ror  (beyond  a  local  smoothing  of  the 
flow)  may  be  unwarranted  in  light  of 
these  resolution  tests. 

HORSESHOE  VORTEX  FLOW 

Representative  results  from  solu¬ 
tions  for  three-dimensional  horseshoe 
vortex  flow  at  zero  and  five  degree  in¬ 
cidence  are  shown  in  Figs.  5a- f.  These 
solutions  converged  in  about  80  time  step 
iterations  and  with  a  14x28x14  grid  Fig.  4  -  Streamwise  Velocity  Profiles 

(a*>5  )  required  about  20  minutes  of  CDC  Near  Separation  for  5:1  Ellipse, 

7600  run  time.  A  14x28x10  grid  was  used  ct*5°,  Re=400. 

for  a*0°.  In  Figs.  5a-b,  vector  plots  of 
velocity  in  a  plane  one  grid  point  away 

from  the  no-slip  endwall  surface  are  shown.  Here,  a  saddle-point  type  of  flow  separa¬ 
tion  is  evident  upstream  of  the  leading  edge  and,  in  the  case  of  a*5°,  toward  the 

(upper)  high  pressure  surface  of  the  ellipse.  The  remaining  plots  in  Fig.  5  show  flow 
velocities  in  the  plane  z=0.15,  located  approximately  in  the  center  of  both  the  horse¬ 
shoe  vortex  flow  region  and  the  approaching  boundary  layer.  Here,  the  velocity  vector 
plots  are  not  much  different  from  the  two-dimensional  flow  region  near  z*l  (not  shown), 
and  the  most  significant  feature  of  the  flow  is  seen  in  the  contours  of  velocity  nor¬ 
mal  to  the  endwall  (and  to  the  page)  in  Figs.  5e-f.  A  strong  downward  flow  toward 
the  endwall  is  present  near  the  leading  edge  (behind  the  saddle-point  separation), 
with  maximum  downward  velocity  of  28  percent  of  the  freestream  reference  velocity 
for  a«0°  and  32  percent  for  a*5  .  Further  results  for  a*5  are  shown  in  Figs.  6a-f. 

2 

Contours  of  total  pressure  loss  coefficient  Ap  /(c  U  /2)  for  z-0.15  are  shown  in 
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Fig.  6a.  In  Fig.  6b,  Che  velocity  in  the  stagnation  plane  normal  to  and  containing 
the  ellipse  leading  edge  again  shows  the  strong  downward  flow  toward  the  endwall  near 
the  leading  edge.  The  cross-flow  velocity  in  outflow  planes  is  shown  in  Figs.  6c-d. 

A  moderately  strong  secondary  flow  pattern  (peak  velocity  20%  of  u  )  indicative  of  a 
streamwise  corner  vortex  is  clearly  in  evidence  near  the  suction  surface  in  Fig.  6d. 
Finally,  a  "limiting"  surface  velocity  vector  plot  and  contours  of  pressure  coeffi¬ 
cient  are  shwon  in  Figs.  6e-f  for  the  surface  one  grid  point  away  from  the  surface  of 
the  ellipse  ("unwrapped"  to  lie  in  a  plane).  Here  there  is  evidence  of  the  distorted 
stagnation  line  near  the  endwall,  flow  toward  the  endwall  near  the  stagnation  line, 
and  finally  the  formation  of  a  streamwise  vortex  visible  mainly  on  the  low  pressure 
or  suction  surface  near  the  endwall. 

CONCLUDING  REMARKS 

Although  no  other  analytical  results  or  experimental  measurements  of  the  three- 
dimensional  horseshoe  vortex  flow  are  available  for  comparison,  the  present  computed 
results  are  consistent  with  flow  visualization  studies  of  related  leading  edge  vortex 
flows.  Computed  results  in  three  dimensions  have  illuminated  several  aspects  of  the 
flow  structure.  For  flows  having  nontrivial  endwall  boundary  layer  development  up¬ 
stream  of  the  leading  edge,  the  horseshoe  vortex  structure  consists  of  an  inertially- 
dominated  rotational  flow  except  very  near  the  solid  boundaries,  where  viscous  effects 
occur  within  thin  layers  generated  locally  by  the  leading  edge  interaction.  The  ro¬ 
tational  inviscid  portion  of  the  vortex  structure  both  upstream  of  the  leading  edge 
and  in  the  downstream  corner  flow  region  scales  with  the  approaching  endwall  boundary 
layer  thickness.  Viscous  effects  are  confined  to  thin  regions  near  the  surfaces, 
having  a  thickness  which  is  apparently  not  strongly  dependent  on  the  approaching 

boundary  layer.  The  overall  flow  consists  of  a  saddle-tvpe  flow  separation  on  the 
endwall  upstream  of  the  leading  edge,  a  strong  spanwise  flow  toward  the  endwall  near 
the  leading  edge,  and  streamwise  vortices  in  the  corner  region  downstream  of  the 
leading  edge.  The  strength  of  the  streamwise  corner  vortex  is  significantly  increased 
on  the  suction  surface  corner  for  flow  at  nonzero  incidence.  Regarding  the  present 
use  of  "zone  embedding"  and  interactive  boundary  conditions  to  minimize  the  computed 
flow  region,  the  general  conclusion  is  that  while  boundary  crnditions  cannot  be 
treated  with  complete  rigor  when  located  within  an  elliptic  region  of  nonuniform 
flow,  careful  treatment  can  lead  to  quite  reasonable  results  which  appear  completely 
adequate  for  the  present  goal  of  horseshoe  vortex  analysis. 
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Fig.  5 


-  Detail  of  Computed  Velocity  for  5:1  Ellipse  with  **0  and  5  ,  Re*400,  M*0.2. 
(a,b):  the  plane  of  grid  points  adjacent  to  endwall  surface;  (c,d):  the 
plane  2*0.15  (near  center  of  approaching  boundary  layer);  (e,f):  contours 
of  velocity  w  normal  to  endwall  in  the  plane  z*0.15. 
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Detail  of  Horseshoe  Vortex  Flow  for  5:1  Ellipse  with  s«5  ,  Re«A00,  M«0.2. 
(a)-Total  pressure  loss  coefficient  in  "0.15  plane;  (b)-Velocitv  in  plane 
containing  leading  edge;  (c,d)-Secondary  velocity  in  outflow  surfaces; 
(e.f)-Velocity  and  contours  of  pressure  coefficient  in  "unwrapped"  surface 
of  points  adjacent  to  ellipse. 


